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Abstract 

We present a preliminary inclusive search for strange pentaquark production in e + e~ interactions 
at a center-of-mass energy of 10.58 GeV using 123 fb^ 1 of data collected with the BABAR detector. 
We look for the states that have been reported previously: the + (154O), interpreted as a ududs 
state; and the S (1860) and H°(1860), candidate dsdsu and uss(uu + dd) states, respectively. 
In addition we search for other members of the antidecuplet and corresponding octet to which 
these states are thought to belong. We find no evidence for the production of such states and 
set preliminary limits on their production cross sections as functions of cm. momentum. The 
corresponding limits on the + (154O) and £ (1860) rates per e + e~ — > qq event are well below 
the rates measured for ordinary baryons of similar mass. 



Submitted to the 32 International Conference on High-Energy Physics, ICHEP 04, 
16 August— 22 August 2004, Beijing, China 



Stanford Linear Accelerator Center, Stanford University, Stanford, CA 94309 



Work supported in part by Department of Energy contract DE-AC03-76SF00515. 



The BABAR Collaboration, 

B. Aubert, R. Barate, D. Boutigny, F. Couderc, J.-M. Gaillard, A. Hichcur, Y. Karyotakis, J. P. Lees, 

V. Tisscrand, A. Zghichc 

Laboratoire de Physique des Particules, F-74941 Annecy-le-Vieux, France 

A. Palano, A. Pompili 
Universitd di Bari, Dipartimento di Fisica and INFN, 1-70126 Bari, Italy 

J. C. Chen, N. D. Qi, G. Rong, P. Wang, Y. S. Zhu 
Institute of High Energy Physics, Beijing 100039, China 

G. Eigcn, I. Oftc, B. Stugu 
University of Bergen, Inst, of Physics, N-5007 Bergen, Norway 

G. S. Abrams, A. W. Borgland, A. B. Breon, D. N. Brown, J. Button-Shafer, R. N. Cahn, E. Charles, 

C. T. Day, M. S. Gill, A. V. Gritsan, Y. Groysman, R. G. Jacobsen, R. W. Kadel, J. Kadyk, L. T. Kcrth, 
Yu. G. Kolomensky, G. Kukartsev, G. Lynch, L. M. Mir, P. J. Oddonc, T. J. Orimoto, M. Pripstcin, 

N. A. Roe, M. T. Ronan, V. G. Shelkov, W. A. Wenzcl 

Lawrence Berkeley National Laboratory and University of California, Berkeley, CA 94720, USA 

M. Barrett, K. E. Ford, T. J. Harrison, A. J. Hart, C. M. Hawkes, S. E. Morgan, A. T. Watson 
University of Birmingham, Birmingham, B15 2TT, United Kingdom 

M. Fritsch, K. Goetzen, T. Held, H. Koch, B. Lewandowski, M. Pelizaeus, M. Steinke 
Ruhr Universitdt Bochum, Institut fur Experimentalphysik 1, D-44^80 Bochum, Germany 

J. T. Boyd, N. Chevalier, W. N. Cottingham, M. P. Kelly, T. E. Latham, F. F. Wilson 
University of Bristol, Bristol BS8 1TL, United Kingdom 

T. Cuhadar-Donszelmann, C. Hearty, N. S. Knccht, T. S. Mattison, J. A. McKenna, D. Thiessen 
University of British Columbia, Vancouver, BC, Canada V6T 1Z1 

A. Khan, P. Kyberd, L. Teodorescu 
Brunei University, Uxbridge, Middlesex UB8 3PH, United Kingdom 

A. E. Blinov, V. E. Blinov, V. P. Druzhinin, V. B. Gohibev, V. N. Ivanchenko, E. A. Kravchenko, 
A. P. Onuchin, S. I. Serednyakov, Yu. I. Skovpen, E. P. Solodov, A. N. Yushkov 

Budker Institute of Nuclear Physics, Novosibirsk 630090, Russia 

D. Best, M. Bruinsma, M. Chao, I. Eschrich, D. Kirkby, A. J. Lankford, M. Mandelkern, R. K. Mommsen, 

W. Roethel, D. P. Stoker 

University of California at Irvine, Irvine, CA 92697, USA 

C. Buchanan, B. L. Hartfiel 

University of California at Los Angeles, Los Angeles, CA 90024, USA 

S. D. Foulkes, J. W. Gary, B. C. Shcn, K. Wang 
University of California at Riverside, Riverside, CA 92521, USA 



2 



D. del Re, H. K. Hadavand, E. J. Hill, D. B. MacFarlane, H. P. Paar, Sh. Rahatlou, V. Sharma 
University of California at San Diego, La Jolla, CA 92093, USA 

J. W. Berryhill, C. Campagnari, B. Dahmes, O. Long, A. Lu, M. A. Mazur, J. D. Richman, W. Verkerke 
University of California at Santa Barbara, Santa Barbara, CA 93106, USA 

T. W. Beck, A. M. Eisner, C. A. Heusch, J. Kroseberg, W. S. Lockman, G. Nesom, T. Schalk, 
B. A. Schumm, A. Seiden, P. Spradlin, D. C. Williams, M. G. Wilson 

University of California at Santa Cruz, Institute for Particle Physics, Santa Cruz, CA 95064, USA 

J. Albert, E. Chen, G. P. Dubois-Felsmann, A. Dvoretskii, D. G. Hitlin, I. Narsky, T. Piatenko, 

F. C. Porter, A. Ryd, A. Samuel, S. Yang 

California Institute of Technology, Pasadena, CA 91125, USA 

S. Jayatilleke, G. Mancinelli, B. T. Meadows, M. D. Sokoloff 

University of Cincinnati, Cincinnati, OH 45221, USA 

T. Abe, F. Blanc, P. Bloom, S. Chen, W. T. Ford, U. Nauenberg, A. Olivas, P. Rankin, J. G. Smith, 

J. Zhang, L. Zhang 

University of Colorado, Boulder, CO 80309, USA 

A. Chen, J. L. Harton, A. Soffcr, W. H. Toki, R. J. Wilson, Q. Zcng 
Colorado State University, Fort Collins, CO 80523, USA 

D. Altenburg, T. Brandt, J. Brose, M. Dickopp, E. Feltresi, A. Hauke, H. M. Lacker, R. Muller-Pfefferkorn, 
R. Nogowski, S. Otto, A. Petzold, J. Schubert, K. R. Schubert, R. Schwierz, B. Spaan, J. E. Sundermann 

Technische Universitat Dresden, Institut fur Kern- und Teilchenphysik, D-01062 Dresden, Germany 

D. Bernard, G. R. Bonncaud, F. Brochard, P. Grenier, S. Schrenk, Ch. Thiebaux, G. Vasileiadis, M. Verderi 

Ecole Poly technique, LLR, F-91128 Palaiseau, France 

D. J. Bard, P. J. Clark, D. Lavin, F. Muheim, S. Playfcr, Y. Xic 
University of Edinburgh, Edinburgh EH9 3JZ, United Kingdom 

M. Andreotti, V. Azzolini, D. Bettoni, C. Bozzi, R. Calabrese, G. Cibinctto, E. Luppi, M. Negrini, 

L. Picmontcse, A. Sarti 

Universitd di Ferrara, Dipartimento di Fisica and INFN, 1-44100 Ferrara, Italy 

E. Treadwell 

Florida A&M University, Tallahassee, FL 32307, USA 

F. Anulli, R. Baldini-Ferroli, A. Calcaterra, R. de Sangro, G. Finocchiaro, P. Patteri, I. M. Peruzzi, 

M. Piccolo, A. Zallo 

Laboratori Nazionali di Frascati dellTNFN, 1-00044 Frascati, Italy 

A. Buzzo, R. Capra, R. Contri, G. Crosetti, M. Lo Vetere, M. Macri, M. R. Monge, S. Passaggio, 

C. Patrignani, E. Robutti, A. Santroni, S. Tosi 

Universitd di Genova, Dipartimento di Fisica and INFN, 1-16146 Genova, Italy 

S. Bailey, G. Brandenburg, K. S. Chaisanguanthum, M. Morii, E. Won 
Harvard University, Cambridge, MA 02138, USA 



3 



R. S. Dubitzky, U. Langenegger 
Universitdt Heidelberg, Physikalisches Institut, Philosophenweg 12, D-69120 Heidelberg, Germany 

W. Bhimji, D. A. Bowerman, P. D. Dauncey, U. Egede, J. R. Gaillard, G. W. Morton, J. A. Nash, 

M. B. Nikolich, G. P. Taylor 

Imperial College London, London, SW7 2AZ, United Kingdom 

M. J. Charles, G. J. Grcnier, U. Mallik 
University of Iowa, Iowa City, IA 52242, USA 

J. Cochran, H. B. Crawley, J. Lamsa, W. T. Meyer, S. Prell, E. I. Rosenberg, A. E. Rubin, J. Yi 
Iowa State University, Ames, IA 50011-3160, USA 

M. Biasini, R. Covarelli, M. Pioppi 
Universita di Perugia, Dipartimento di Fisica and INFN, 1-06100 Perugia, Italy 

M. Davier, X. Giroux, G. Grosdidicr, A. Hocker, S. Laplace, F. Le Diberder, V. Lcpclticr, A. M. Lutz, 
T. C. Petersen, S. Plaszczynski, M. H. Schunc, L. Tantot, G. Wormser 

Laboratoire de VAccelerateur Lineaire, F-91898 Orsay, France 

C. H. Cheng, D. J. Langc, M. C. Simani, D. M. Wright 
Lawrence Livermore National Laboratory, Livermore, CA 94550, USA 

A. J. Bevan, C. A. Chavez, J. P. Coleman, I. J. Forster, J. R. Fry, E. Gabathuler, R. Garnet, 
D. E. Hutchcroft, R. J. Parry, D. J. Payne, R. J. Sloane, C. Touramanis 

University of Liverpool, Liverpool L69 72E, United Kingdom 

J. J. Back, 1 C. M. Cormack, P. F. Harrison, 1 F. Di Lodovico, G. B. Mohanty 1 
Queen Mary, University of London, El 4NS, United Kingdom 

C. L. Brown, G. Cowan, R. L. Flack, H. U. Flaecher, M. G. Green, P. S. Jackson, T. R. McMahon, 

S. Ricciardi, F. Salvatore, M. A. Winter 

University of London, Royal Holloway and Bedford New College, Egham, Surrey TW20 OEX, 

United Kingdom 

D. Brown, C. L. Davis 
University of Louisville, Louisville, KY 40292, USA 

J. Allison, N. R. Barlow, R. J. Barlow, P. A. Hart, M. C. Hodgkinson, G. D. Lafferty, A. J. Lyon, 

J. C. Williams 

University of Manchester, Manchester M13 9PL, United Kingdom 

A. Farbin, W. D. Hulsbergen, A. Jawahery, D. Kovalskyi, C. K. Lae, V. Lillard, D. A. Roberts 
University of Maryland, College Park, MD 20142, USA 

G. Blaylock, C. Dallapiccola, K. T. Flood, S. S. Hertzbach, R. Roller, V. B. Koptchev, T. B. Moore, 

S. Saremi, H. Staengle, S. Willocq 

University of Massachusetts, Amherst, MA 01003, USA 

x Now at Department of Physics, University of Warwick, Coventry, United Kingdom 



4 



R. Cowan, G. Sciolla, S. J. Sekula, F. Taylor, R. K. Yamamoto 
Massachusetts Institute of Technology, Laboratory for Nuclear Science, Cambridge, MA 02139, USA 

D. J. J. Mangeol, P. M. Patcl, S. H. Robertson 
McGill University, Montreal, QC, Canada H3A 2T8 

A. Lazzaro, V. Lombardo, F. Palombo 
Universita di Milano, Dipartimento di Fisica and INFN, 1-20133 Milano, Italy 

J. M. Bauer, L. Cremaldi, V. Eschenburg, R. Godang, R. Krocgcr, J. Rcidy, D. A. Sanders, D. J. Summers, 

H. W. Zhao 

University of Mississippi, University, MS 38611, USA 

S. Brunet, D. Cote, P. Taras 
Universite de Montreal, Laboratoire Rene J. A. Levesque, Montreal, QC, Canada H3C 3J1 

H. Nicholson 

Mount Holyoke College, South Hadley, MA 01015, USA 

N. Cavallo, 2 F. Fabozzi, 2 C. Gatto, L. Lista, D. Monorchio, P. Paolucci, D. Piccolo, C. Sciacca 
Universita di Napoli Federico II, Dipartimento di Scienze Fisiche and INFN, 1-80126, Napoli, Italy 

M. Baak, H. Bultcn, G. Raven, H. L. Snoek, L. Wildcn 

NIKHEF, National Institute for Nuclear Physics and High Energy Physics, NL-1009 DB Amsterdam, 

The Netherlands 

C. P. Jessop, J. M. LoSecco 
University of Notre Dame, Notre Dame, IN 46556, USA 

T. Allmcndinger, K. K. Gan, K. Honscheid, D. Hufnagcl, H. Kagan, R. Kass, T. Pulliam, A. M. Rahimi, 

R. Ter-Antonyan, Q. K. Wong 

Ohio State University, Columbus, OH 43210, USA 

J. Brau, R. Frey, O. Igonkina, C. T. Potter, N. B. Sincv, D. Strom, E. Torrence 
University of Oregon, Eugene, OR 91403, USA 

F. Colecchia, A. Dorigo, F. Galeazzi, M. Margoni, M. Morandin, M. Posocco, M. Rotondo, F. Simonetto, 

R. Stroili, G. Tiozzo, C. Voci 

Universita di Padova, Dipartimento di Fisica and INFN, 1-35131 Padova, Italy 

M. Benayoun, H. Briand, J. Chauveau, P. David, Ch. de la Vaissiere, L. Del Buono, O. Hamon, 
M. J. J. John, Ph. Leruste, J. Malcles, J. Ocariz, M. Pivk, L. Roos, S. T'Jampens, G. Therin 

Universites Paris VI et VII, Laboratoire de Physique Nucleaire et de Hautes Energies, F-15252 Paris, 

France 

P. F. Manfredi, V. Re 
Universita di Pavia, Dipartimento di Elettronica and INFN, 1-21100 Pavia, Italy 
2 Also with Universita della Basilicata, Potenza, Italy 



5 



P. K. Bchera, L. Gladncy, Q. H. Guo, J. Panetta 
University of Pennsylvania, Philadelphia, PA 19104, USA 



C. Angelini, G. Batignani, S. Bcttarini, M. Bondioli, F. Bucci, G. Calderini, M. Carpinelli, F. Forti, 
M. A. Giorgi, A. Lusiani, G. Marchiori, F. Martinez- Vidal, 3 M. Morganti, N. Ncri, E. Paoloni, M. Rama, 

G. Rizzo, F. Sandrelli, J. Walsh 

Universitd di Pisa, Dipartimento di Fisica, Scuola Normale Superiore and INFN, 1-56127 Pisa, Italy 

M. Hairc, D. Judd, K. Paick, D. E. Wagoner 
Prairie View A&M University, Prairie View, TX 77446, USA 

N. Danielson, P. Elmer, Y. P. Lau, C. Lu, V. Miftakov, J. Olsen, A. J. S. Smith, A. V. Tclnov 
Princeton University, Princeton, NJ 08544, USA 

F. Bellini, G. Cavoto, 4 R. Faccini, F. Ferrarotto, F. Ferroni, M. Gaspero, L. Li Gioi, M. A. Mazzoni, 
S. Morganti, M. Pierini, G. Piredda, F. Safai Tehrani, C. Vocna 

Universitd di Roma La Sapienza, Dipartimento di Fisica and INFN, 1-00185 Roma, Italy 

S. Christ, G. Wagner, R. Waldi 
Universitdt Rostock, D- 18051 Rostock, Germany 

T. Adye, N. De Groot, B. Franck, N. I. Geddes, G. P. Gopal, E. O. Olaiya 
Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, 0X11 OQX, United Kingdom 

R. Aleksan, S. Emery, A. Gaidot, S. F. Ganzhur, P.-F. Giraud, G. Hamel de Monchenault, W. Kozanecki, 
M. Legendre, G. W. London, B. Mayer, G. Schott, G. Vasseur, Ch. Ycche, M. Zito 

DSM/Dapnia, CEA/Saclay, F- 9 1191 Gif-sur-Yvette, France 

M. V. Purohit, A. W. Weidemann, J. R. Wilson, F. X. Yumiceva 
University of South Carolina, Columbia, SC 29208, USA 

D. Aston, R. Bartoldus, N. Berger, A. M. Boyarski, O. L. Buchmueller, R. Claus, M. R. Convery, 
M. Cristinziani, G. De Nardo, D. Dong, J. Dorfan, D. Dujmic, W. Dunwoodie, E. E. Elsen, S. Fan, 
R. C. Field, T. Glanzman, S. J. Gowdy, T. Hadig, V. Halyo, C. Hast, T. Hryn'ova, W. R. Innes, 
M. H. Kelsey, P. Kim, M. L. Kocian, D. W. G. S. Lcith, J. Libby, S. Luitz, V. Luth, H. L. Lynch, 
H. Marsiske, R. Messner, D. R. Muller, C. P. O'Grady, V. E. Ozcan, A. Pcrazzo, M. Perl, S. Petrak, 
B. N. Ratcliff, A. Roodman, A. A. Salnikov, R. H. Schindler, J. Schwicning, G. Simi, A. Snyder, A. Soha, 
J. Stelzer, D. Su, M. K. Sullivan, J. Va'vra, S. R. Wagner, M. Weaver, A. J. R. Weinstein, 
W. J. Wisniewski, M. Wittgen, D. H. Wright, A. K. Yarritu, C. C. Young 

Stanford Linear Accelerator Center, Stanford, CA 94309, USA 

P. R. Burchat, A. J. Edwards, T. I. Meyer, B. A. Petersen, C. Roat 
Stanford University, Stanford, CA 94305-4060, USA 

S. Ahmed, M. S. Alam, J. A. Ernst, M. A. Saeed, M. Saleem, F. R. Wappler 

State University of New York, Albany, NY 12222, USA 

3 Also with IFIC, Instituto de Fisica Corpuscular, CSIC-Universidad de Valencia, Valencia, Spain 
4 Also with Princeton University, Princeton, USA 



6 



W. Bugg, M. Krishnamurthy, S. M. Spanier 
University of Tennessee, Knoxville, TN 37996, USA 

R. Eckmann, H. Kim, J. L. Ritchie, A. Satpathy, R. F. Schwitters 
University of Texas at Austin, Austin, TX 78712, USA 

J. M. Izen, I. Kitayama, X. C. Lou, S. Ye 

University of Texas at Dallas, Richardson, TX 75083, USA 

F. Bianchi, M. Bona, F. Gallo, D. Gamba 
Universitd di Torino, Dipartimento di Fisica Sperimentale and INFN, 1-10125 Torino, Italy 

L. Bosisio, C. Cartaro, F. Cossutti, G. Delia Ricca, S. Dittongo, S. Grancagnolo, L. Lanceri, P. Poropat, 5 

L. Vitalc, G. Vuagnin 

Universitd di Trieste, Dipartimento di Fisica and INFN, 1-34127 Trieste, Italy 

R. S. Panvini 

Vanderbilt University, Nashville, TN 37235, USA 

Sw. Banerjee, C. M. Brown, D. Fortin, P. D. Jackson, R. Kowalewski, J. M. Roney, R. J. Sobie 
University of Victoria, Victoria, BC, Canada V8W 3P6 

H. R. Band, B. Cheng, S. Dasu, M. Datta, A. M. Eichcnbaum, M. Graham, J. J. Hollar, J. R. Johnson, 
P. E. Kutter, H. Li, R. Liu, A. Mihalyi, A. K. Mohapatra, Y. Pan, R. Prepost, P. Tan, J. H. von 

Wimmersperg-Toeller, J. Wu, S. L. Wu, Z. Yu 

University of Wisconsin, Madison, WI 53706, USA 

M. G. Greene, H. Neal 
Yale University, New Haven, CT 06511, USA 



5 Deceased 



7 



1 Introduction 



Recently several experimental groups have reported observations of a new, manifestly exotic (B=l, 
S=l) baryon resonance, called the + (154O) [l]-[7], with an unusually narrow width (r <8 MeV/c 2 ) 
for a particle in this mass range that has open channels to strong decay. Also, the NA49 experiment 
reported evidence for an additional narrow exotic (B=l, Q=S= —2) state, called S , as well as 
corresponding S° state, with masses close to 1862 MeV/c 2 [8]. More recently the HI collaboration 
reported a narrow (T < 30 MeV/c 2 ) exotic charmed (B=l, C= —1) resonance, 0°, with a mass 
of 3099 MeV/c 2 . The simplest quark assignments consistent with the quantum numbers of , 
S and 0° are (ududs), (dsdsu) and (ududc), respectively; therefore these observed states are 
regarded as pentaquark candidates. 

These results have prompted a surge of pentaquark searches in experimental data of many 
kinds, mostly with negative results [10]. Several theoretical models [11, 12, 13] have been proposed 
to describe possible pentaquark structure. They predict that the lowest mass states containing u, d 
and s quarks should occupy a spin- 1/2 antidecuplet and octet, as illustrated in Fig. 1. Predictions 
for the masses of the unobserved N and U states vary; Mjy might be anywhere between ~150 
MeV/c 2 below Me and ~ Ms, and Mjj anywhere between ~ M@ and ~150 MeV/c 2 above Ms- 
In order to distinguish pentaquarks from ordinary baryons, we adopt a notation that has names 
corresponding to ordinary baryons with similar s quark content, plus a subscript 5. For example, 
the 3$ pentaquark (dss(uu + dd)) corresponds to the ordinary 3~ (dss) baryon. 




Figure 1: Quark structure of the antidecuplet (annuli) and octet (filled circles) that are generally 
assumed for the lowest-mass pentaquarks. Strangeness increases in the vertical direction and the 
third component of isospin in the horizontal. 
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Although experiments with a baryon in the beam or the target might have some advantage in 
pentaquark production, e + e~ interactions are also known for democratic production of hadrons. 
Baryons with nonzero beauty, charm, strangeness (up to three units) and/or orbital angular mo- 
mentum have been observed with production rates that appear to depend on the mass and spin, 
but not quark content. If pentaquarks are produced similarly, then one might expect a pentaquark 
rate as high as that for an ordinary baryon of the same mass and spin, i.e. about 8 x 10 -4 0§ and 
about 4 x E§ per e + e~ — ► hadrons event at y/s = 10.58 GeV [15]. Decays of B mesons are 
also known to produce a rather high rate of baryons and so provide another fertile hunting ground. 
In both cases we expect any set of pentaquark states forming a multiplet of approximately equal 
mass to be produced at roughly equal rates, so that access to all narrow states is available. 

Here we describe a program of inclusive searches for the reported states and also the other 
members of the hypothetical antidecuplet and octet in data from the BABAR experiment, which 
include both e + e~ — > qq, q = udsc, events and the production of T(4S) mesons, which decay into 
BB pairs. We have concentrated on decay modes involving strange particles and protons, which 
are easily reconstructed in our detector. All states except the A% have nonzero strangeness, and 
the N$ and U$ have hidden strangeness divided in some way between the octet and antidecuplet 
states. Since the production mechanism, and hence the momentum spectra, for pentaquarks is 
unknown, we consider differential production cross sections da /dp* per unit momentum p* in the 
e + e~ center-of-mass (cm.) frame, which are to first order independent of any model, and applicable 
to the sum of all production processes at yfs = 10.58 GeV. After describing some common aspects 
of the analyses in section 2, we describe searches in three classes of decay modes in the next three 
sections. Upper limits on pentaquark production are discussed in section 6, and we summarize the 
results in section 7. 

2 Data and Simulation 

The data sample used in the analysis comprises the runs of the BABAR experiment through 2003, 
and amounts to an integrated luminosity of 123 fb _1 , roughly 90% (10%) of which was recorded on 
(below) the T(4S) resonance at a cm. energy of y/s = 10.58 (10.54) GeV. The BABAR detector 
is described in detail in Ref. [14]. Here we use charged tracks reconstructed in the 5-layer silicon 
vertex detector (SVT) and the 40-layer drift chamber (DCH), and identified using a combination 
of energy loss measured in these two subdetectors and Cherenkov angles measured in the detector 
of internally reflected Cherenkov radiation (DIRC). The charged particle momentum resolution is 
given by (5pt/pt) 2 = (0.0013pr) 2 + 0.0045 2 , where pr is the momentum transverse to the beam 
axis in GeV/ c. Loose particle identification is required on most particles, giving in each case nearly 
100% efficiency with a sizable background reduction. Tighter requirements are made where noted, 
giving identification efficiencies ranging from 80-99% and background rejection factors of 12-100. 
Photons are detected by a CsI(Tl) electromagnetic calorimeter (EMC) with an energy resolution 
oia{E)/E = 0.023 ■ {E /GeV)' 1 ^ 0.019. Our invariant mass resolution for a given decay mode is 
better than that of any experiment reporting the observation of a pentaquark candidate with that 
mode available. 

Minimal hadronic event selection is performed, as we wish to be as inclusive as possible and 
maintain maximum signal efficiency. We simply require three reconstructed charged tracks in the 
event. The vast majority of the events are hadronic; lepton-pair events are expected to contribute 
negligible background in the signal regions; two-photon events and events with a hard initial-state 
photon contribute only at low momenta. The acceptance for a pentaquark signal from any of these 
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sources is well understood; if a signal is found, we can attempt to isolate the source with cuts on 
the event properties. 

Samples of simulated e + e~ — > qq, where q = udsc, are generated using the JETSET [16] Monte 
Carlo generator combined with a detailed simulation of the BABAR detector, in order to study 
backgrounds and "control" particles - known resonances with similar masses and decay modes 
to the pentaquark in question. By comparing with the same particles reconstructed in the data, 
we verify the simulation of the invariant mass resolutions and biases of the detector, an essential 
ingredient in the evaluation of any new signal, or limit where no signal is seen. 

In order to evaluate pentaquark reconstruction efficiencies, we simulate pentaquark signals 
within our standard software with special parameter settings of the JETSET generator. For each 
pentaquark state we take an existing baryon as a stand-in, change its mass and width to match 
the desired pentaquark properties, and force a particular decay mode, with all other JETSET 
parameters left unchanged. The generator thus gives the stand-in particles a momentum spectrum 
based on the pentaquark mass, but the stand-in flavor. We make no attempt here to conserve flavor 
in the production process, as we are only concerned with providing an appropriate number and 
distribution of additional particles in the event. We do not know the production mechanism for 
such five-quark states; however, to make the measurement we only need the efficiency as a function 
of the relevant variables, in this case momentum and polar angle in the laboratory frame, so any 
smooth distribution that covers the entire range of these variables is sufficient. Effects such as the 
type, multiplicity and proximity of other particles in the event can be studied by using different 
stand-in baryons; such effects are found to be smaller than the other systematic uncertainties. The 
stand-in particles used to generate our signal samples are summarized in Table 1. 

Table 1: Baryons used in the simulation to generate pentaquark signal samples. 



Pentaquark 


Existing 


Mass 


Width 


#of 


Signal mode 


Baryon Used 


(MeV/c 2 ) 


(MeV/c 2 ) 


Events 




A+ 


1540 




120,000 




A++ 


1862 




120,000 






1862 
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3 Search for the 0+(154O) -> pK° s 

We reconstruct 65" candidates in the pK° s decay mode, where K° s — > 7r~7r + , with selection 
criteria designed for high efficiency and low bias against any production mechanism. A 
sample of K® candidates is obtained from all pairs of oppositely charged tracks passing 
loose pion identification requirements that pass within 6 mm of each other. Each pair is 
required to: have a total momentum vector extrapolating within 6 mm (32 mm) of the 
interaction point (IP) in the plane transverse to (along) the beam direction; have a positive 
flight distance, defined as the projection on the candidate momentum direction of a vector 
from its point of closest approach to the beam axis to its decay point; and have an invariant 
mass within 10 MeV/c 2 of the nominal K° s mass. We also require the helicity angle 9 H of the 
candidate, defined as the angle between the ir + and the reconstructed it + it~ flight direction 
in the 7r + 7r~ rest frame, to satisfy \cos6h\ < 0.8, eliminating background from A decays 
and photon conversions and reducing combinatoric background. 

We combine the surviving K° s candidates with identified p and p tracks that extrapolate 
within 15 mm (10 cm) of the IP in the plane transverse to (along) the beam direction. The 
simulated signal reconstruction efficiency varies with p*, the pK® momentum in the e + e~ 
cm. frame, from 13% at low p* to 22% at high p*. The invariant mass distribution of these 
pairs is shown in Fig. 2. There is a clear peak at 2285 MeV/c 2 from A£ — > pK® but no 
other sharp structure. The A+ peak (upper right plot) shows a mass resolution of better 
than 6 MeV/c 2 and contains roughly 52,000 entries, demonstrating our sensitivity to the 
presence of a narrow resonance. The lower plots zoom in on the mass ranges 1400-1600 and 
1600-1800 MeV/c 2 . The 9^ has been reported at values between 1520 and 1540 MeV/c 2 . 
There is no enhancement in our data anywhere in this region. The pK® decay mode is also 
possible for the E$ states, whose masses might be expected to be between that of the 65 
and about 2000 MeV/c 2 depending on their strange quark content. There is no sign of a 
narrow resonance anywhere in this region. 

We consider several additional criteria that might enhance a pentaquark signal. More 
stringent requirements on the flight distance of the K° s candidate give a much cleaner K® 
sample at the expense of efficiency at low momentum. This enhances the A£ signal to 
noise, but does not reveal any additional structure. If a pentaquark is produced in an e + e~ 
annihilation event, then there must also be an antibaryon (or anti-pentaquark) in the event, 
among whose decay products is either an antiproton or antineutron. In the case of the Q% 
decaying to pK° S) the K° s must have been a K° rather than a K°, and there must be a 
compensating particle in the event with strangeness — 1, which might often be a K~ . We 
consider all subsets of the pK° s combinations for which there is also at least one loosely 
identified K~ and/or p track in the event. In all cases the data quantity is reduced, the 
A+ signal is still visible and there is no sign of a pentaquark peak. Following work by the 
CLAS and NA49 experiments [3, 8], we reconstruct a possible iV*(2400) — > K~K®p using 
the selected pK® pairs and an additional loosely identified K~ candidate. We observe no 
iV* peak in our data. Considering various K~ K° s p mass ranges, the pK® mass distribution 
for K~KgP masses near 2400 MeV/c 2 shows no pentaquark signal, and is indistinguishable 
from the pK® mass distributions in nearby K~K®p mass ranges. Requiring in addition a 
recoil antiproton yields similar results. 
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Figure 2: Distribution of the pK® invariant mass for combinations satisfying all the criteria de- 
scribed in the text. The same data are plotted four times in different pK® mass regions. 
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In each of the above cases we split the data into ten bins of p*, in order to enhance our 
sensitivity to any production mechanism that gives a momentum spectrum different from 
that of our background. The bins are 500 MeV/ c 2 wide and cover the p* range from zero to 
5 GeV/c 2 , the kinematic limit for a particle of mass 1700 MeV/c 2 . The background is much 
smaller at higher p*, so we are more sensitive to mechanisms that produce harder spectra. 
In no case is there any sign of a pentaquark signal. 

We quantify these null results for the Gg~ , assuming a mass of 1540 MeV/c 2 . In order 
to reduce model dependence, we consider the ten p* bins noted above and fit a signal plus 
background function to the pK® invariant mass distribution in each bin. The natural width 
of the 65" has not been measured; the best upper limit, T < 8 MeV/c 2 , is larger than our 
detector resolution, so we must consider the range of widths up to this value. We use a 
P-wave Breit-Wigner lineshape multiplied by a phase-space factor and convolved with a 
resolution function derived from the A£ data and simulation. The latter is a sum of two 
Gaussian distributions with a common center and an overall rms ranging from 2.5 MeV/c 2 
at low p* to 1.8 MeV/c 2 at high p*; this is narrower than for the A£ due to the proximity 
of the 0^(1540) to threshold. For the natural width we consider two possiblities, 1 MeV/c 2 
and 8 MeV/c 2 , corresponding to a very narrow state and the upper limit, respectively. The 
fits are performed over a wide range, from threshold to 1800 MeV/c 2 , and the background 
function is a seventh-order polynomial times a threshold factor. 

For the case of our nominal selection criteria (Fig. 2), we obtain the signal yields shown 
in Fig. 3. In all p* bins the fit quality is good and the signal is consistent with zero. There 
is no positive trend in the data, and the roughly symmetric scatter of the points about zero 
indicates low momentum- dependent bias in the background function. We consider systematic 
effects in the fitting procedure by varying the signal and background functions and fit range; 
changes in the fitted signal are negligible compared with the statistical uncertainties. Other 
selection criteria yield similar results. Since the nominal selection gives the smallest absolute 
uncertainties after efficiency corrections, we use it to set upper limits on the production 
cross section in section 6. Varying the mass assumed for the Bg" has effects consistent with 
statistics; in no case do we observe a signal, and the uncertainties shown in Fig. 3 are typical. 
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Figure 3: The Q§ signal yields extracted from the fits to the pK$ invariant mass distributions, 
assuming a 9^" mass of 1540 MeV/c 2 and width T = 1 MeV/c 2 (left) and T = 8 MeV/c 2 (right). 
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4 Search for S 5 — >S 7r and Hg — > 3 7r + 

We next search for pentaquark states decaying into a £~ and a charged pion, where S~ — > 
yl 7r _ and yi° — > p7r~, including the reported S 5 (1862) and £;5(1862). We first reconstruct 
A — > p7r candidates from all pairs of charged tracks satisfying loose proton and pion 
requirements. Efficient and unbiased selection criteria are again applied: the tracks must 
pass within 6 mm of each other; the candidate have a positive flight distance from the IP; and 
it must have an invariant mass within 10 MeV/c 2 of the nominal A mass. These candidates 
are combined with an additional negatively charged track passing loose pion identification 
requirements to form E~ candidates. These candidates are required to form a good vertex, to 
have a positive flight distance from the IP, and to have an invariant mass within 20 MeV/ c 2 
of the nominal S~ mass. Furthermore the flight distance of the A candidate from the A°7r~ 
vertex must be positive. Finally, we combine the S~ candidates with an additional charged 
track consistent with coming from the IP and passing loose pion identification requirements. 
The cosine of the angle between the reconstructed E~ trajectory, extrapolated back to the 
IP, and the additional track is required to be less than 0.998. This last cut is especially 
important, since the S~ is charged and has a long lifetime; if it has a long flight distance, 
it can produce a reconstructed track that, if combined with itself, forms a false peak in the 
invariant mass distribution. The simulated signal reconstruction efficiency varies from 6.5% 
at low p* to 12% at high p*. 

The invariant mass distributions for E~tt~ and E~tt + are shown in Figs. 4 and 5, respec- 
tively. In Fig. 5 there are clear peaks as expected for the £r*°(1530) and £; ( ?(2470) baryons, 
but no other structure is visible. In Fig. 4 there is no visible sharp structure at all. 




Figure 4: E tt invariant mass distribution. 
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Figure 5: S~ir + invariant mass distribution. 

As in the preceding section, we divide the S 5 candidates into ten bins of p* and find no 
sign of a pentaquark signal in any bin. We fit a similar signal plus background function to 
the S~tt~ invariant mass distribution in each bin. The resolution function is derived in this 
case from the S*°(1530) and .5 , °(2470) data and simulation, and is described by a Gaussian 
with an rms of 8 MeV/c 2 . For the Breit-Wigner width we consider two possiblities, 1 and 
18 MeV/c 2 , corresponding to a very narrow state and the upper limit on the S 5 width, 
respectively. The fits are performed over the mass range from 1760 to 1960 MeV/c 2 , and the 
background function is a seventh-order polynomial. 

Fixing the S 5 mass to 1862 MeV/c 2 , we obtain the signal yields shown in Fig. 6. In all 
p* bins the fit quality is good across the full mass range and the signal is consistent with zero. 
Systematic uncertainties on the fitting procedure are again found to be negligible compared 
with the statistical uncertainties, and variations of the S 5 mass and selection criteria give 
consistent results. 

5 Search for S~, SJ, i\T° and AT+ in the A°K and U°K modes 

We search for pentaquark states decaying into a A°K or E°K final state, where K is a 
charged or neutral kaon, E° — > A 0r y and A — > pn~ . These final states give us access to the 
Sf, S®, N® and N$ states (see Fig. 1). The has been reported with a mass of 1862 
MeV/c 2 , and we expect a quite similar mass for the Sf; since our decay modes include two 
strange particles, we are sensitive only to N® and N§ states with substantial ss content, and 
these might be expected somewhere between 1500 and 1900 MeV/c 2 . 

We first reconstruct A — > pn~ candidates from all pairs of charged tracks in which one 
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Figure 6: The S§ signal yields extracted from the fits to the S~ir~ invariant mass distributions, 
assuming a E§ mass of 1862 MeV/c 2 and width T = 1 (dashed) and T = 18 MeV/c 2 (solid). 

satisfies tight proton identification requirements, and the other loose pion identification. The 
pairs are required to form a good vertex, to have an angle between their flight direction (line 
from the IP to the vertex) and total momentum at the vertex less than 200 mrad, and 
to have an invariant mass within 20 MeV/c 2 of the nominal A mass. To reconstruct S° 
candidates, The A candidates are combined with neutral (not associated with any charged 
track) calorimetric energy deposits of at least 80 MeV that do not pair with any other neutral 
deposit to form a ir° candidate. We require a mass difference M pnl — M p7T within 20 MeV/c 2 
of the nominal M s o — M^o value. Similarly, K° s candidates are reconstructed from pairs 
of oppositely charged tracks forming a good vertex, having an angle between their flight 
direction and total momentum less than 30 mrad, and having an invariant mass within 4.5 
MeV/c 2 of the nominal K® mass. Charged kaon candidates are required to be consistent 
with coming from the IP and to pass tight kaon identification requirements. The simulated 
signal reconstruction efficiencies vary from 0.5% for the £°K® mode to 10% for the A°K ± 
modes at low p*, and increase to 3% and 25%, respectively, at high p* . 

We find a substantial p* dependence in the structure of the invariant mass distributions, 
and so we show them in four p* bins, for A°K + , A°K~ and A°K° S combinations in Figs. 7, 
8 and 9, respectively, and for S°K + , E°K~ and E K® combinations in Figs. 10, 11 and 12, 
respectively. In Fig. 7 there is a peak as expected from yl+(2285) and also some structure in 
the 1950-2150 MeV/c 2 region that can be attributed to A^ decays with one or two missing 
pions. No other narrow peaks are evident. 

In Fig. 8 the i?~ peak is prominent and there is a hint of a broad structure corresponding 
to the £;~(1820), but no other narrow structure is visible. In Fig. 9 we observe the S , °(2472) 
at high momenta, and there are hints of the i; (1820), but no other narrow structure. The 
E° K modes have much lower statistics, but show structure consistent with the corresponding 
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Figure 7: Invariant mass distributions of A°K + for p* in four different ranges. 



A°K mode. There is a hint of the (1950) at low p* in the £°K~ mode. In no case is 
any unexpected narrow structure seen, and all evident structure is consistent with that due 
to known resonances. We have also examined the data in finer p* bins and find no sign of a 
pentaquark signal in any p* range. 

We perform fits in ten p* bins to the A°K~ and A G K Q S mass distributions, in which 
we might expect signals near 1860 MeV/c 2 from the and states, respectively. The 
resolution functions derived from the data and simulation are described by sums of two 
Gaussians with total rms of about 4 MeV/c 2 that depend slightly on p*. For each Breit- 
Wigner width we consider both 1 and 18 MeV/ c 2 , as for the S 5 . The fits are performed over 
the mass range from 1772 MeV/c 2 to 1972 MeV/c 2 in order to exclude known resonances. 
The background function is a first-order polynomial. 

We fix the masses to 1862 MeV/c 2 , plus a shift obtained from the simulation in each 
p* bin that does not exceed 2 MeV/c 2 , and obtain the signal yields shown in Fig. 13. In 
all p* bins the fit quality is good across the full mass range and the signal is consistent 
with zero. Systematic uncertainties on the fitting procedure are again found to be negligible 
compared with the statistical uncertainties, and variations of the mass and selection criteria 
give consistent results. 
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Figure 8: Invariant mass distributions of A°K for p* in four different ranges. 
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Figure 9: Invariant mass distributions of A Kg for p* in four different ranges. 
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Figure 10: Invariant mass distributions of E°K + for p* in four different ranges. 
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Figure 11: Invariant mass distributions of £°K~ for p* in four different ranges. 
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Figure 12: Invariant mass distributions of U°Kg for p* in four different ranges. 
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6 Upper Limits 



For the states reported by previous experiments, there exist specific masses at which to 
search and experimental upper limits on the decay widths. We can therefore calculate upper 
limits on their differential production cross sections, with some assumptions on the branching 
fraction to the mode in which the search was made. We present such limits for the 0^~, , 
~z and 

For the G^~(1540) — > pK° s we take the yields from the fits shown above (see Fig. 3) and 
convert them into cross sections by dividing by the efficiency (including the K s — > tt~tt + 
branching fraction), the integrated luminosity and the p* bin width. The reconstruction 
efficiency is calculated from the simulation and corrected using data, and varies smoothly 
from 13% at low p* to 22% at high p*. It is checked by measuring the differential cross 
section for A£ — > pK® in the combination of qq and T(4S) events represented in our data, 
which is found to be consistent with the appropriate combination of previous measurements. 
If the G>5~ is a ududs pentaquark state, we expect only two possible decay modes, nK + and 
pK°, with very similar Q values and hence roughly equal branching fractions. Assuming 
half of the K° appear as Kg, we arrive at a branching fraction B(0§ — > pK®) = 1/4, and 
we divide by this value to obtain the total ©5" differential production cross section. 

We then derive a conservative upper limit on this cross section in each bin. We consider 
only the physically allowed region, and scale the limit by a factor of (1 + Se/e), where 
5e/e = 0.049 is the sum in quadrature of the relative systematic uncertainties on the efficiency 
and luminosity. The relative luminosity uncertainty is 1%, and that on the efficiency is 
dominated by the 3.2% uncertainty on the reconstruction of pairs of tracks from a displaced 
vertex such as a Kg. The total is nearly independent of p*. 
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Figure 14: The 95% C.L. upper limit on the production cross-section for the 0^ assuming a mass 
of 1540 MeV/c 2 and a natural width T = 1 MeV/c 2 (left) or T = 8 MeV/c 2 (right), as a function 
of cm. momentum p*. 



These upper limits are shown in Fig. 14 and tabulated in Table 2. Two sets of limits 
are shown; one corresponds to a very narrow state, and the other to a width at the current 
experimental upper limit of 8 MeV/c 2 . The limits correspond to the mass of 1540 MeV/c 2 
used in the fits; repeating the analysis at several nearby mass values gives consistent limits. 
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The units are fb per GeV/c, and apply to the sum of all possible production processes. To 
isolate continuum and T(4S) production we divide by the respective cross sections, and the 
corresponding limits on the numbers of pentaquarks per event are given in Table 2. 

Similarly we convert the measured yields described above for the S 5 (1862) — > E~tt~ 
(Fig. 6), S 5 "(1862) -> A°K- (Fig. 13) and H 5 °(1862) -> A°K° S (Fig. 13) decays into cross 
sections. Here we assume a mass of 1862 MeV/ c 2 and present limits for both a very narrow 
hypothesis and for T = 18 MeV/c 2 . The S 5 (1862) — > E~7t~ reconstruction efficiency varies 
smoothly from 6.5% at low p* to 12% at high p*, and has been checked using the observed 
£;*°(1530) — > E~tt + signal. Its relative systematic uncertainty again varies slightly with p*, 
and its average value of 6e/e = 0.064 is larger than for the pK® mode, as there are two 
displaced vertices and more tracks in the decay. The £;^~(1862) — > A°K~ efficiency varies 
from 10% at low p* to 27% at high p*, is checked with the Q~ — > A°K~ signal, and the 
average Se/e = 0.052 is similar to the pK° s mode. The E®{1&Q2) — > A Q K Q S efficiency varies 
from 3.5% at low p* to 13% at high p*, and 5e/e = 0.072 is similar to the S~ir~ mode. 

We assume a S b — ► E~-k~ branching fraction of one-half; for the S§ and we do not 
assume branching fractions into the measured modes, but present limits on the product of 
the production cross section and branching fraction. Again the results are independent of 
the assumed mass; they are shown in Figs. 15 and 16, and tabulated in Tables 3, 4 and 5. 
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Figure 15: The 95% C.L. upper limit on the production cross-section for S 5 assuming a natural 
width of r = 1 (dashed) and T = 18 MeV/c 2 (solid), as a function of cm. momentum p*. 

In order to quote limits on the total production cross section (times branching fraction) 
we must either know the momentum spectrum or believe that it does not vary rapidly on the 
scale of our bin size. Since the latter is true in the simulation, we assign conservative, model- 
independent upper limits on the total numbers of pentaquarks produced per e + e~ — > qq event 
(T(45) decay) by summing each differential cross section over the p* range from zero to the 
kinematic limit for 5.3 GeV jets (B meson decays), taking into account the fact that most 
of the systematic error is common to all bins. Limits are derived from these sums by the 
same method as in each bin, and are listed in Tables 2-5. Any postulated spectrum can be 
folded with our differential limit to obtain a limit on the total cross section assuming that 
spectrum, which will be lower than the model independent limits given in the tables. 
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Figure 16: The 95% C.L. upper limits on the production cross-section for S 5 times its branching 
fraction into A°K~ (left), and times its branching fraction into A Kg (right), assuming a natural 
width of r = 1 MeV/c 2 (solid) and T = 18 MeV/c 2 (dashed), as a function of cm. momentum p*. 



7 Summary 

We have performed a preliminary high statistics search for the reported states (1540), 
S 5 (1862), Sf (1862) and (1862) in e + e~ annihilations, and also for most of the other 
members of the pentaquark octet and anti-decuplet to which they are postulated to belong. 
In all cases we observe clear signals for known baryon resonances that demonstrate sensitivity 
to any new narrow resonances, with invariant mass resolution better than the reported upper 
limits on the widths of the respective states. We find no evidence for the production of such 
states in 123fb _1 of BABAR data. For the reported states, we see no excess at the measured 
invariant mass values and use the reported limits on their widths to set upper limits on the 
their inclusive production in our search modes in both e + e" — > qq events events and T(4S) 
decays. These limits are at the level of 10 -4 -10~ 5 per event, depending on the width assumed, 
and are valid for any narrow state in the vicinity of 1540 or 1860 MeV/c 2 . The searches for 
other members of the multiplets show no evidence for an unknown narrow resonance at any 
mass between threshold for the decay mode used and the kinematic limit. 

In order to limit the production rates we must know the branching fraction of each state 
into the decay mode used in the search. In the case of the G^~, we take this to be 1/4, 
as noted above. For E 5 a similar argument that E~tt~ and E~ K~ dominate and have 
roughly equal branching fractions leads to a branching fraction of 1/2. Taking the upper 
limit widths, we calculate 95% C.L. upper limits on the total production rates of 1.1 x 10 -4 
G5" and 1.0 x 10~ 5 E 5 per e + e~ — > qq event (preliminary); these are roughly a factor of 
eight and four below the typical values measured for ordinary octet and decuplet baryons of 
the same masses of 8 x 10~ 4 and 4 x 10~ 5 , respectively, as illustrated in Fig. 17. 

The situation is more complex for the E§ and E® as the mixing between the members of 
the anti-decuplet and the octet is unknown. In these two cases the branching fraction of an 
antidecuplet state of mass ~1860 MeV/c 2 to the measured mode could be anywhere from 
zero to ~l/3. The best limit that could be set, assuming the value of 1/3, would therefore be 
roughly 15 x 10 -5 Sf or per e + e~ — > qq event, which is well above the "expected" value 
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Table 2: The measured signal yield in each p* bin, assuming a mass of 1540 MeV/c 2 and two 
values of the natural width T in MeV/c 2 . The corresponding upper limits at the 95% C.L. on the 
05~ differential production cross section, the number of produced per e + e~ — ► qq event and the 
number per T(4S) decay. 
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of 4 x 10 5 . The study of additional modes is needed to elucidate the production properties 
of these states. 
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Table 3: The measured S 5 signal yield in each p* bin, assuming a mass of 1862 MeV/c 2 and two 
values of the natural width in MeV/c 2 . The corresponding upper limits at the 95% C.L. on the 
S 5 differential production cross section, the number of S~ 5 produced per e + e~ — ► qq event and 
the number per T(45) decay. 
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Table 4: The measured ST" signal yield in each p* bin, assuming a mass of 1862 MeV/c 2 and two 
values of the natural width in MeV/c 2 . The corresponding upper limits at the 95% C.L. on the E§ 
differential production cross section times its branching fraction B into A°K~, Bxthe number of 
E$ produced per e + e~ — > qq event and .Bxthe number per T(4S) decay. 
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Table 5: The measured signal yield in each p* bin, assuming a mass of 1862 MeV/c 2 and two 
values of the natural width in MeV/c 2 . The corresponding upper limits at the 95% C.L. on the 
differential production cross section times its branching fraction B into A°K®, i?xthe number of 
£5 produced per e + e~ — ► qq event and .Bxthe number per Y(4S) decay. 
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Figure 17: Compilation of baryon production rates in e + e~ annihilation [15] from experiments at 
the Z° (circles) and ~ 10 GeV (squares) as a function of baryon mass. The vertical scale 
accounts for the number of spin and particle+antiparticle states, and the lines are chosen to guide 
the eye. The arrows indicate our preliminary upper limits on spin-1/2 @f and E§ pentaquark 
states, assuming the branching fractions shown, and are seen to lie below the solid line. 



27 



References 

[1] LEPS Collaboration, T. Nakano et ai, Phys. Rev. Lett. 91, 012002 (2003). 

[2] SAPHIR Collaboration, J. Barth et ai, Phys. Lett. B 572, 127 (2003). 

[3] CLAS Collaboration, V. Kubarovsky et al, Phys. Rev. Lett. 92, 032001 (2004). Erra- 
tum; ibid, 049902. 

[4] DIANA Collaboration, V.V. Barmin et al, Phys. Atom. Nucl. 66, 1715 (2003). 

[5] SVD Collaboration, A. Aleev et al, hep-ex/0401024 (2004). 

[6] HERMES Collaboration, A. Airapetian et al, Phys. Lett. B 585, 213 (2004). 

[7] COSY-TOF Collaboration, M. Abdel-Bary et ai, Phys. Lett. B 595, 127 (2004). 

[8] NA49 Collaboration, C. Alt et ai, Phys. Rev. Lett. 92, 042003 (2004). 

[9] HI Collaboration, A. Aktas et ai, hep-ex/0403017 (2004). 
[10] See, e.g., J. Pochodzalla, hep-ex/0406077 (2004) and references therein. 
[11] D. Diakonov, V. Petrov and M.V. Polyakov, Z. Phys. A 359, 305 (1997). 
[12] M. Karliner and H. Lipkin, hep-ph/0402260 (2004). 
[13] R. Jaffe and F. Wilczek, Phys. Rev. Lett. 91, 232003 (2003). 

[14] The BABAR Collaboration, B. Aubert et al, Nucl. Instrum. Methods A 479, 1 (2002). 
[15] Particle Data Group, K. Hagiwara et al, Phys. Rev. D 66, 010001 (2002). 
[16] T. Sjostrand, Comput. Phys. Commun. 82, 74 (1994). 



28 



